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1 Introduction 

An antibody, also known as an immunoglobulin, is a protein used by the immune system to identify, neutralize, 
or kill foreign objects like bacteria, viruses, or pollen which are termed as antigen. The production of antibodies 
is the main function of the immune system. An antigen, when introduced into the body, triggers the production 
of an antibody by immune system which will then kill or neutralize the antigen that is recognized as a foreign 
invader. 

The bio-medical application of antibodies against an effect of toxins associated with various biological threats 
(epidemic outbreaks or bio-terrorism) is well-documented (see, e.g., [1-3]). 

For a long time the main target of antibody design has been the antibody affinity. With progress in bio- 
engineering, many antibodies with different affinity parameters have been generated. However, according to 
Skvortsov and Gray [4] affinity is not a good predictor of protective or therapeutic potential of an antibody. 
In fact, the treatment effect of an antibody can be described by a parameter which includes the reaction rates 
of the receptor-toxin-antibody (RTA) kinetics and relative concentration of reacting species. As a result, any 
given value of this parameter determines a range of antibody kinetic properties and its relative concentration 
in order to achieve a desirable therapeutic effect. 

The model considered by Skvortsov and Gray is a model of a well-mixed solution of toxin, antibody, and cells 
and neglects diffusion fluxes of interacting species. Diffusion fluxes are significant especially when the process 
of RTA interaction is limited by diffusion. Skakauskas et al. [5] examined numerically a RTA interaction model 
taking into account diffusion of all species in the case where a spherical cell is embedded into an initially 
uniformly distributed toxin-antibody solution which occupies a large volume (compartment) lying between the 
cell and external surface. Initial values of species and their values on the external surface were assumed to be 
the same for all times. In this case fluxes of toxin, antibody and their complex across the external surface are 
not zero. Some numerical results of the evaluation of an antibody treatment efficiency parameter are given in 
this paper. 

In the present paper by using the same model we study the influence of RTA kinetic parameters and diffusivity 
of toxin, antibody, and their complex on the behavior of the antibody protection parameter and concentrations 
of species in more detail. 

The paper is organized as follows. In Section 2 we introduce the reaction-diffusion model for RTA interaction. 
Numerical results are presented in Section 3. Summarizing remarks given in Section 4 conclude the paper. 



2 The model 

We study a case of a spherical cell embedded into a toxin-antibody solution which occupies an extracellular 
domain SI lying between the cell and an external surface and use notations of paper [5]: 
p - spherical radius, 

S c = {p : p = p c } - the surface of the spherical cell, p c is its radius, 

S e = {p : p = p e } - the surface of the external sphere (external surface of SI), p e is its radius, 

SI = {p : p G (p c , p e )} - the extracellular domain, 

ro — the concentration of receptors on the cell surface, 

0(t,p) - the fraction of the toxin-bound receptors, 

rgO - the concentration of the toxin-bound receptors (confined to S c ) , 

r (l — 9) - the concentration of the free receptors, 

ut, ua, and uc - the concentrations of toxin, antibody, and toxin-antibody complex, respectively, 
«p u^, Uq - the initial concentrations, 
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Kt, ka, and kq - the diffusivity of the toxin, antibody, and toxin-antibody complex, respectively, 

fei, k—i - the forward and reverse constants of the toxin-antibody reaction rate, 

k 2 and k-2 - the forward and reverse constants of the toxin and receptor binding rate, 

&3 - the rate constant of the toxin internalization, 

d n - the outward normal derivative on S e or S c , 

d t =d/dt, 

A = p 2 -^-(p 2 Tr) ~ the Laplace operator, 
dp op 

ip(t) - the antibody protection factor (a relative reduction of toxin inside a cell due to application of 
antibody) . 

Dynamics of the concentrations ut, ua, uc, and 9 can be described by the following equations: 



' d t UT — —kiUrUA + k-\Uc + ktAut, p € fl, t > 0, 
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(3) 



d t uc = kiu T u A - k-iuc + kcAmc, 

U C =0, p = p e , t > 0, 
d n uc = 0, p = p c , t > 0, 

,«c|t=o = 0, pGfi. 



peU, t > 0, 



(4) 



The initial and boundary conditions for the system above correspond to a case where initially the toxin and 
antibody are distributed uniformly in the extracellular domain Q. Values of all species on the outer boundary 
of 12 for all times and their initial values are assumed to be the same. In particular, zero value of the toxin- 
antibody complex is used for initial time and for all times on the outer boundary of H We stress that in this 
case the fluxes of all species are not zero on the outer boundary S e of Q. 

Eqs. (l)-(4) can be presented in non-dimensional form by using scales of r* (time), I (length), and u„ 
(concentration). By substituting variables x = Ix, t = r»f, r = lu*f n , ut — u„ut, ua — u*ua, uc — u*uc, 



uto 
K A = 



= lt»«y, UAQ = U*U A , kl = T*U»fcl, k 2 = T*U*k 2 , fc_l = T„fc_l, fc_ 2 = T*fc_ 2 , fc 3 = T„fc 3 , K T = T*K T l 2 , 

t*kaI~ 2 , kc = t*kcI~ 2 into (l)-(4) we can deduce the same system, but only in the non-dimensional 



variables. Therefore, for simplicity in what follows, we treat system (l)-(4) as non-dimensional. 

The main antibody treatment efficiency parameter is the antibody protection factor (a relative reduction of 
toxin attached to a cell due to application of antibody) which can be defined by the following expression [4,5]: 



Isc d \u'X>a dS 
f Sc °\u° A =o dS' 



(5) 



By definition < tp < 1. The lower the value of tp the more profound is therapeutic effect of antibody 
treatment. 



3 Numerical results 

We treated system (l)-(4) numerically for the spherically symmetric domain, p e (p c ,p e ), and t > with an 
implicit finite-difference scheme. Our selection of the values of parameters was motivated by the values available 
in the literature [3,5-7] with the extended range to allow exploration and illustration of the various transport 
and kinetics regimes that are possible in the RTA system. We employ the following data that were used in the 
most calculations in [5,8]: w* = 6.02 • 10 13 cm' 3 , r* = 1 s, r = 1.6 • 10 4 / S c , where 1.6 • 10 4 is the total number of 



receptors of the cell, Z = 10 2 cm, S c = 47rp 2 = 47r • 10 6 cm 2 , r = 2.115 • 10 3 . The standard non-dimensional 
values of the other parameters are the following: 

'fci = 1.3 • 10~ 2 , k-i = 1.4- 10~ 4 , 
k 2 = 1.25 • 10" 2 , fc_ 2 = 5.2- 10- 4 , fc 3 = 3.3-10- 5 , 
< k t = 10~ 2 , k a = 10~ 2 , K C = 10~ 2 , (6) 

Pc = 10- 1 , /5 e = 2, 

u o =1; u o =05 

These values correspond to the ricin and 2B11 mono-clonal antibody interaction. If values of k\, k 2 , ka, 
kc, and kt differ from those given in (6), they are specified in the legends of plots. 

As we indicated in the Introduction, the main purpose of our study was to estimate the effect of diffusive 
and kinetic parameters of species on the behavior of concentrations of species and protective properties of an 
antibody against a toxin. Results of numerical solving of system (l)-(4) are presented in Figs. 1-7. 
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Fig. 1. Influence of the external radius p e = 2 (solid line) and 5 (dashed line) and the toxin diffusivity kt ■ W 2 (1), 
5 • 1(F 3 (2), 1CP 3 (3), 1CF 4 (4) on the cell protection characteristic, ip, in the case of u° T = 0.6. 





The plots of tp in Fig. 1 depict the dependence of the antibody protection factor on the radius p e of the 
external surface S e and toxin diffusivity kt- Parameter tp increases with kt growing, but its behavior for large 
values of kt is non-monotonic. For large values of kt, parameter tp grows as p e decreases. But for small values 
of kt its behavior is different. For example, if kt = 10~ 3 , then values of ip for p e = 5 are larger than those for 
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Fig. 3. Effect of the toxin diffusivity kt : 10~ 2 (solid line), 5 • 10 -3 (dashed line), 10~ 3 (dash-dotted line) and 
parameter fei : 1.3 • 1CT 2 (1), 2 x 1.3 • 1CT 2 (2), 4 x 1.3 • 10" 2 (3) on the cell protection function ip. 
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Fig. 4. Effect of the toxin diffusivity k t : KT 2 (solid line), 5 • KT 3 (dashed line), 2.5 • 1CT 3 (dash-dotted line) and 
parameter k 2 : 1.25 • KT 2 (1), 2 x 1.25 • 1CT 2 (2), 4 x 1.25 • 1CT 2 (3) on the cell protection function ip. 



p e = 2 if t < 1400 s approximately. But if kt < 10~ 4 , then, for all t, values of ip for p e = 5 are larger than 
those for p e = 2 (see curves 3 and 4). 

Fig. 2 illustrates the dependence of ip on the diffusivity ka of the antibody. The curves in this figure depict 
the increase of ip as ka decreases and non-monotonic time evolution of ip for small values of ka- Moreover, in 
the case of small antibody diffusivity, ka = 10~ 4 , values of ip for p e = 2 are larger than those for p e = 5. But 
in the case of large antibody diffusivity, ka = 10 _1 , values of ip for p e = 2 are smaller than those for p e = 5. 
only if t < 1000 s. For t > 1000 s they behave vica versa. 

Figs. 3 and 4 exhibit the dependence of ip on diffusivity Kt, forward constant k\ of the toxin and antibody 
reaction rate, and forward constant ki of the toxin and receptor binding rate, respectively. Fig. 3 demonstrates 
the decrease of ip as k\ increases. But different values of k\ do not change the monotonic behavior of all curves 
in time. From Fig. 4 we see the non-monotonic behavior of ip as &2 increases. Moreover, ip increases with &2 
increasing. The bottom of the hollow in Fig. 3 is located lower than that in Fig. 4. One can see in Fig. 4 
that the effect of toxin diffusivity variation on protection factor is sensitive to changes of parameter ki . Let us 
compare the minimal values of protection factor. In the case of k^ = 1-25 • 10~ 2 , the minimum of ip is about 
0.72 at kt = 10~ 2 , 0.62 at kt — 5 • 10~ 3 and 0.5 at kt — 2.5 • 10~ 3 , while the corresponding values of ip are 
about 0.86, 0.84 and 0.836 in the case of &2 = 5 • 10~ 2 (curves 1 and 3). 

Numerical experiments show that diffusivity kq practically does not influence the time evaluation of ip. 

The plots of ut in Fig. 5 depict the dependence of the toxin concentration at the cell surface on the diffusivity 
kt and radius p e of the external surface S e . For any value of p e , ut decreases with kt decreasing. For large 



values of Kt, function UT(t,p c ) grows as p e decreases. But, for small values of kt, its behavior is different. For 
example, for kt = 10~ 3 , values of UT(t,p c ) for p e = 5 are larger than those for p e = 2 only if t < 500 s (see 
curves 3). Our calculations show that influence of Kc on the behavior of ux{t, p c ) is insignificant. We observed 
the non-monotonic behavior of UT{t,p c ) for small kc, but difference between its steady-state value and value 
at the bottom of the hollow is very small (of order 10~ 3 ). 

Calculations show that uc(t,p c ) grows with kc decreasing. The behavior uc{t, p c ) is monotonic for kc < 
5 • 10~ 2 . Its values are smaller than initial ones of toxin for Kp e [5 ■ 10~ 3 , 5 • 10~ 2 ]. But uc(t, p c ) can reach a 
relatively large steady-state value for small kc while steady-state values of ut and ua are smaller than their 
initial values. For example, the steady-state value of uq on S c is equal to 2.2 for kc = 10~ 3 , ka = kt = 10~ 2 . 
Derivatives of ut and ua with respect to p on S e arc of order 0.3 while derivative of uc on S e is of order —3. 
This means that uc increases faster towards the cell than ut and ua decays in the same direction. 

Curves in Fig. 6 depict the dependence of 6 on kt and p e for u A = 1 and u A = 0. In the case where 
the antibody is absent values of 9 practically do not depend on diffusivity kt (see the bullets marked curve). 
9 decreases with kt decreasing. For any p e , function 9 grows as kt increases. If Kt € 5 • [10 -3 , 10~ 2 ], then 
values of 9 for p e = 2 are larger than those for p e = 5. But for small values of kt its behavior is different. For 
example, if kt = 10 -3 , values of 9 for p e — 5 are larger than those for p e = 2 only for about t < 1300 s. This 
behavior is similar to those of ut and tp. 

Two curves in Fig. 7 illustrate the non-monotonic behavior of derivatives duA(t, p e )/dp and duc{t, p e )/dp 
for small toxin diffusivity (kt = 10~ 3 ). For kt = 10~ 2 their behave is monotonic. 

4 Concluding remarks 

To conclude the paper we summarize results of study. The receptor-toxin-antibody interaction is studied 
numerically by using a model proposed in [5]. The model includes "bulk" reaction of toxin and antibody, 
surface binding of toxin and cell receptors, and diffusion of all species. The main results of the numerical study 
are the following: 

1. The evolution of concentrations of some species (toxin and toxin-bound receptors) and of the antibody 
protection factor for some cases (large toxin diffusivity, small antibody diffusivity, and large forward constant 
of the toxin-receptor binding rate) is non-monotonic 

2. The influence of small or large values of Kt, ka, and on the behavior of UT(t,p c ), 9(t) and ip(t) is 
profoundly different in the cases of small or large p e . 

3. The effect of kc on the evolution of ut, 9, and ip was found to be insignificant. 
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